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CHAPTER 1 



INTRODUCTION 
1 . 1 GENERAL 

One of the most important characteristics of a nuclear 
reactor is the effect of temperature on reactivity. As the 
temperature of the moderator, fuel, or reflector are 
individually changed, the effect on the reactivity may be 
quite different. Because it is difficult to mock up a 
temperature profile in a low power reactor, it is common to 
discuss temperature effects in terms of the isothermal 
temperature coefficient of reactivity where isothermal implies 
that all parts of the reactor are at a constant and uniform 
temperature. Comparison of the measured and calculated iso- 
thermal temperature coefficient aids progress toward under- 
standing the physical processes taking place in the reactor. 

This dissertation is concerned with two aspects of 
thermal effects on reactivity. The first is an experimental 
study of the isothermal temperature coefficient of reactivity 
of the Cornell University Zero Power Reactor, a light-water 
moderated, 2.1% enriched, UO 2 fueled, aluminum clad reactor 
having nominal water- to-fuel ratios of 1:1, 1.5:1, 2:1, 3:1, 
and 4:1.—/ (See Appendix A for fuel and lattice parameters.) 

— / S. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores," CURL-7, p. 220 Cornell 
University Reactor Laboratory (1964) . 
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A theoretical prediction is given in Chapter 2, and the experi- 
mental techniques and results are discussed in Chapter 3. 

The second aspect of this work singled out the temperature 
dependence of the disadvantage factor as one contribution to 
the temperature dependent behavior of the reactor. Chapter 4 
discusses the theoretical basis of the use of foils to measure 
disadvantage factors. Experimental methods and results are 

2 / 

discussed in Chapter 5 and compared with predictions by MacVean.— 
Chapter 6 summarizes the results of both aspects of this 
work and suggests further experimental studies. 

1.2 ISOTHERMAL TEMPERATURE COEFFICIENT OF REACTIVITY 
The isothermal temperature coefficient a is defined as 
dp/d9 where p is excess reactivity and 9 is the temperature 
of all parts of the core. Isothermal coefficients of reactiv- 
ity have previously been measured for cores similar to the 
Cornell critical assembly. In these studies, temperature 
dependent excess reactivities were determined either by 
measuring periods for a given control rod configuration or 
by measuring the change in critical position for a calibrated 
control rod. 

2 / 

— C. R. MacVean, "A Simplified Cell Theory Applied to the 
Calculation of Thermal Neutron Spectra in Light Water 
Lattices," CURL-10, p. 145 Cornell University Reactor 
Laboratory (1964) . 
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The Westinghouse temperature coefficient experiment s -^- 1 — ' 
were made on aluminum clad, 1 . 3%, enriched metal and oxide 
fuels of both 0.600 and 0 .387-inch diameters. Reactivity 
values were determined by measurement of periods longer than 
100 seconds. Water- to- fuel ratios from 1.5:1 to 5:1 were 
studied over a temperature range of 20-50°C. 

The Yankee cores, also studied by Westinghouse,— ^ had 
2.2:1 and 2„9:1 water- to-fuel ratios for 2.7% enriched, 
stainless steel clad, 0.300-inch diameter U0 9 fuel. Reactiv- 
ity determination was by period measurement with various 
control rod configurations. The temperature range covered 
was 22-75°C. 

The Cornell University Zero Power Reactor was studied 
over a temperature range 5-70°C for water- to-fuel ratios of 
1:1, 2:1, and 3:1. After clearly establishing that the 4:1 

— J. R. Brown, et al^ , "Reactor Properties of Water Moderated 

Slightly Enriched Uranium Lattices," WAPD-117, p. 117 
Westinghouse Atomic Power Division (1954) . 

— J. R. Brown, e_t al^., "Kinetic and Buckling Measurements on 

Lattices of Slightly Enriched Uranium or UO 2 Rods in 
Light Water," WAPD-176, p. 39 Westinghouse Atomic Power 
Division (1958) . 

— P. W. Davison, et al., "Yankee Critical Experiments-Measure- 

ments on Lattices of Stainless Steel Clad Slightly Enriched 
Uranium Dioxide Fuel Rods in Light Water " YAEC-94, n„ "" 
Yankee Atomic Electric Company (1959) . 
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lattice had a positive isothermal coefficient of reactivity 
with no indication of going negative, the study of this core 
was terminated at a maximum temperature of 45°C. The 
accuracy of the values of a for both the Westinghouse and 
Cornell experiments is about + 0.03 4/°C, and this accuracy 
is nearly independent of the value of a. 

13 TEMPERATURE DEPENDENT DISADVANTAGE FACTORS 

The neutron density disadvantage factor 6^ and the 

experimentally measured disadvantage factor & eX p are defined 

by the ratios 6 = t^/N^ and 5 = C™ /CT^ where N and (T are 

J n exp 

respectively the spatially averaged neutron density and foil 
count rate per unit volume; superscripts m and f indicate 
respectively the moderator and fuel regions of the reactor. 

Disadvantage factor measurements have previously been 
made on cores similar to the Cornell critical assembly. The 
Brookhaven measurements ^ 1 — 1 — 1 were made with dysprosium 

— H. Kouts, "Intra-cell Flux Traverses," BNL-1783, 

Brookhaven National Laboratory (1953) . 

— ■ H. Kouts, "Intracell Flux Traverses and Thermal Utilization 
for 1.15% Enriched Uranium Rods in Ordinary Water," 
BNL-1982, Brookhaven National Laboratory (1954). 

— ! Go A. Price, "Thermal Utilization Measurement," BNL-1992, 
Brookhaven National Laboratory (1954) . 
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pellets distributed through the fuel and moderator regions. 

The pellets were 1/16-inch diameter, 0.007-inch thick foils 
made from a mixture of Dy^O^ with aluminum or plastic. 

Numerical integration of the detailed flux maps gave average 
count rates for each region which were converted to experi- 
mental disadvantage factors. An uncertainty of about Tk was 
assigned from data scatter; statistical uncertainty in 
counting was about 17> . 

The Westinghouse disadvantage factor measurements -— 

were made with uranium sector foils cut and placed in the 

fuel and moderator regions of the reactor. The foil was 

235 

highly enriched U dispersed in aluminum to match the fuel 
235 

U atom density. These foils were gamma-counted for 
fission product activity. The experimental disadvantage 
factor was obtained directly as the ratio of count rates per 
unit volume of the moderator and fuel sector foils. The 
uncertainty in the measurement was about 37>. 

— Ho J. Kouts, et al., "Intracell Flux Traverses and Thermal 

Utilization, 1.0277 o Enriched Uranium Rods in Light Water," 
BNL-1796, Brookhaven National Laboratory (1954) . 

— ^ H. C. Honeck, "The Calculation of Thermal Utilization and 
Disadvantage Factors in Uranium/Water Lattices," Nuc . 

Sci. Eng. L8, 49 (1964). 
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No previously published works indicate that the tempera- 
ture dependence of the disadvantage factor has been studied 

experimentally. MacVean's recently developed method of cal- 

13/ 

culating thermal lattice parameters — gives the temperature 
dependence of the disadvantage factor. Since the calculated 
changes are of the order of 10 /°C an experimental method 

was required which was more accurate than those previously 
used . 

Experimental disadvantage factors were measured for 
three water- to-fuel ratios over a temperature range 5-65°C. 

An integral foil counting method was used employing a 
0.600 0.001-inch diameter foil in the fuel and an appro- 

priately shaped moderator foil. Five-mil thick foils of 
dysprosium (5.24% dysprosium in aluminum), manganese (P-metal), 
and vanadium were chosen because their cross sections are 
nearly 1/v. Two counting methods were used. One was to 
dissolve the manganese foils in nitric acid and gamma-count 
aliquots of the solutions in a well-type scintillation 
counter. The other method was to beta-count the solid foils 
in a proportional flow counter. Results of the two methods 
agreed within the statistical errors which were about 0.7%. 

— ! A. Z. Franz, "Measurements of Thermal Utilization Resonance 
Escape Probability, and Fast Fission Factor of Water 
Moderated Slightly Enriched Uranium Lattices," WAPD-134, 
Westinghouse Atomic Power Division (1955) . 
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— — ' A. Z. Kranz and G. G. Smith, "A Second Report on Measure- 
ments of f, p, and e of Water Moderated Slightly En- 
riched Uranium Lattices," WAPD-151, Westinghouse Atomic 
Power Division (1956) . 

— ' C. R. MacVean, footnote 2, pp . 143-145. 



CHAPTER 2 



THEORETICAL CAL CULATION OF THE I SOTHE RMAL T EMPE R ATU RE 
CO EFFICIENT OF RE ACT IV ITY 

2.1 INTRODUCTION 

A calculation of the isothermal temperature coefficient 
of reactivity was made using two group, two region pertur- 
bation theory. The presentation here follows closely that 
14 / 

of Larrimore.—— 

Section 2.2 presents and discusses the two group, two 
region neutron balance equations and the associated boundary 
conditions necessary for the solution. The bilinearly 
weighted, first order reactivity equations from perturbation 
theory are given for the core, reflector, and interface 
regions of the reactor. In these latter equations, tempera 
ture derivatives of the reactor parameters are weighted by 
the fluxes, adjoint fluxes, flux gradients, and adjoint flux 
gradients. The temperature derivatives are discussed in 
Sec. 2.3. Utilization of these equations for machine cal- 
culations is discussed in Sec. 2.4. 

— ^ J. A. Larrimore, "Temperature Coefficients of Reactivity 
in Homogenized Thermal Nuclear Reactors," Ph.D. Thesis, 
pp . 185-211, M.I.T. (September, 1962). 
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2 . 2 WO GROUP, WO REGION EQUATIONS 
2,2.1 St eady State Neutron Bal an ce Equati on s 

The steady state neutron balance equations for the core 
region of the reactor are: 



^ D 1c V2-: 1c ) cp lc +VG2 fis Cp 2c 


= 0 


(2.1) 


p Tc^lc + ( ^ D 2 c 72 ' 2 2c ) 9 2c 


= 0 . 


(2.2) 



Analogous equations for the reflector are 

(D lr y 2 ~z lr ) cp lr = 0 (2.3; 

^lr^lr ^2r^ ” 2 2r^ ^ 2 t = ^ *(2.4) 

where: subscript 1 refers to the fast group, subscript 2 

refers to the thermal group, subscript c denotes the core, 

subscript r denotes the reflector, D is the diffusion constant 

for a given group and region, 2^ is the fast group removal 

cross section, Z £ . is the fission cross section, 2 0 is the 

thermal absorption cross section, ) is the average number 

of neutrons per thermal fission, e is the fast fission factor, 
2 

^ is the Laplacian differential operator, and p is the 
resonance escape probability. 

The adjoint equations for the core and reflector are 
given by 

< D lc^- Z lc> + Phc^c ' 0 < 2 - 5 > 

u ' 5Z fis <P lc + < D 2 c v2 -"2c> t 2c * 



0 



( 2 . 6 ) 
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(2.7) 




0 



( 2 . 8 ) 



where the superscript + indicates the adjoint flux functions 

Boundary conditions at the outer extrapolated boundary 

require that the fast and slow fluxes and adjoint fluxes 

vanish, At the core-reflector interface the boundary con 

ditions require that fast and slow fluxes and adjoint fluxes 

match (i.e cp^ = c P^ r > c Pq^ c = 9^) and similar equations with 

the 1 replaced by 2), and that the fast and slow currents 

and adjoint currents match (i.e. D- cp. = D. cd. where i 
J v ic ic lr lr 

is either 1 or 2, and the same expressions are written for 
the the adjoint fluxes). 

22,2 Bilinearly Weighted Reactivity Relations 

Once the fast and slow fluxes and adjoint fluxes for 
the core and reflector are known, they can be used in the 
standard two group, two region, first-order perturbation 
theory reactivity equations. Perturbations of the parameters 
are equated to temperature derivatives of the parameters, 
and the reactivity effects are calculated. 

The following expression is given by Larrimore-^-^ for 
the core contribution to the isothermal temperature co- 
efficient, a : 

’ c 



15 / 



J A Larrimore, footnote 14, p. 191. 
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a = 
c 



/ d - 



*lc v€ -fis tp 2cl 



/ 



dr [-W> lc • 



dDn , dD~ , dZ-, 

”35“ Vnc “V^c ' ”35“ V<P 2 c " ^lc “30“ ‘Pic 



+ <f>t„ V€ 



dZ.. 

f IS 



lc Vfe “HT” Cp 2c ^ 1+1 2 c ~3© 



9 



+ d( P Z lc> 



cp 



lc 



dZ 0 

Cp 2c “S© 2 cp 2c 



1 



(2.9) 



where dir is the volume element, the subscript c on the 

integral sign implies integration over the core of the 

reactor, and 0 is temperature. 

Similarly, the reflector contribution, a^ s is given 
16/ 

by Larrimore-^-i' as 

Jdr <p{ c vcS £is <p 2 \ 1 | Jdr [-V< r • 

c / ' r 



a = 
r 



c 

dD 



dDr 



dS 



35” ^^lr "V c P2r * ”30”^/ cp 2r ~ 9 lr “3©” ^lr 



+ 9 2r “3?T ^lr " ^r “35^ ^20 



dZ 9 
+ 2r 



( 2 . 10 ) 



-j r / 

— J. A. Larrimore, footnote 14, p. 194. 



Sec . 2.2 



12 



Finally, Larrimore gives for the interface effect,™^ 




+ cp 2r cp lc^ p2 lc" rj lr^ ) 

+ 9 t r cp 2c veZ f . s - 4 r 9 2c (S 2c -2 2r )l . 

-J mt 



[dV/de] 



( 2 . 11 ) 



where all the functions in the second set of brackets are 
evaluated at the core-reflector interface, and dv/d0 is the 
change of core volume with temperature. For an infinite 
cylinder, dV/dQ is given by 

dV/d0 = 2irr 2 P a] _/3 (2.12) 



where r is the core radius, and is the volumetric 

efficient of thermal expansion of aluminum. 



co- 
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J. A. Larrimore, footnote 14, p. 209. 
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The total isothermal temperature coefficient for the 
reactor is the sum of the core, reflector, and interface 
contributions : 



a = a +a +a. . (2.13) 

c r l 

2.3 TEMPERATURE DERIVATIVES OF THE NUCLEAR PARAMETERS 
In order to calculate the various expressions for the 
isothermal temperature coefficient, it is necessary to have 
relations for the temperature derivatives of the various 
nuclear parameters, viz., v , e , p , and the various D's and 
Z's. 

2.3.1 Reflector Group Constants 

The temperature derivative of the macroscopic cross 
section Z is easily found by differentiating 

2 = N a o , (2.14) 

where is the number density of atoms whose microscopic 
cross section is a to give 

(l/2)dZ/d0 = (l/N )dN / d0 + (l/a)dc/d0 . (2.15 ; 

d 3. 

Since N is proportional to the density d of the material, 

3. 

it follows that 



(l/N a )dN a /d0 = -p , 



(2.16) 
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where [3 is the volumetric coefficient of thermal expansion 
for the material. By substitution in Eq, (2.15), the tempera- 
ture derivative of 2 becomes 

(l/Z)dZ/d© = - (3+ (. 1 / a) do/d0 . (2.17) 

For the fast group, it was assumed that the microscopic 
cross section does not change with temperature. This 
assumption leaves only a density effect term for the fast- 
group macroscopic cross section temperature derivative. For 

the reflector we therefore have 

(l/S lr ) d ^ir/d0 ' -e H2 o . (2.18) 

The slow microscopic cross section for water o 2r is 

- 1/2 

proportional to 1/v; hence it is proportional to 9 . 

Therefore, 

(l/a 2r )do 2r /d0 = -0,5/9 (2,19) 

where 9 is in °K. For the slow group in the reflector, 

Eq. (2.17) therefore becomes 

(l/2 2r )d2 2r /d9 = “P H q-0.5/9 - (2.20 

The diffusion constant D can be written 

D = l/(3i, tr ) , (2.21) 

where is the appropriate transport cross section. 
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Logarithmic differentiation gives 

(l/D)dD/d© = +p-da tr /d© . (2.22) 

18 / 

Larrimore—— ' gives the following relationship for the 
logarithmic slope of the thermal transport cross section 
of water: 

(l/°tr) da tr/ d0 = ' n//Q * (2.23) 

where n is given for 20, 50, and 100°C as 0.473, 0.471, and 
0.469 respectively. Using the value for 50 o C, Eq. (2.22) 
becomes 

(l/D 2r )dD 2r /d© = ^0,471/0 . (2.24) 

Assuming again that the microscopic cross section of the 
fast group does not change with temperature, Eq. (2.22) 
for the fast group becomes 

< 1/D lr> dD lr /d ® " \ 0 (2 ' 25 > 

2.3.2 Core Group Constants 

In the core region of the reactor, the fast group 
macroscopic cross section 2 c can be written 

E c = (V f +V m ) -1 (V f Z f +V n '^ m ) , (2.26) 

18 / 

— — ' J. A. Larrimore, footnote 14, pp 83-84 
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where and V m are the volumes of the fuel and moderator 
region of the cell, 2^ is the fuel cross section, and 2 m 
is the moderator cross section. 

For the fast group, the microscopic cross sections are 
assumed not to vary with temperature, and the expansion of 
the fuel is taken to be negligible compared to the expansion 
of the cladding. The logarithmic temperature derivative of 
2-^ c then becomes 

(l/Zi c )dZi c /d0 - V m (V f +V m )- 1 (S lr /2 lc ) ^2 B a 1 /3)-6 h J , 



(2.27) 

where (3^1 is the volumetric coefficient of thermal expansion 
of aluminum, the cladding material. The last term in 
Eq . (2.27) arises because the density effect on the moderator 
is the combined expansion effects of the water and aluminum 
For the reactor cell, it is only the water cross-sectional 
area change per unit length of the reactor cell which con- 
tributes to the water density change. This area change is 
due to the expansion of aluminum structural members. The 
2 ^ai/ 3 term accounts for this area change per unit length. 
Using Eq. (2.21) where, for the fast group, 

2 tr = (V f 4V m )" 1 (V m 2 1 ^ r +V f 2^ r ) , (2.28) 



it can be shown that the logarithmic temperature derivative 
of D lc is 
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(l/D lc )dD lc /d0 = ( D i c /Di r )(V m )(V f +V ni )" 1 ^ H 0 “23 a 1 /3) . 

(2.29) 



The temperature derivatives for the slow group core 

parameters were taken from calculations by MacVean using the 

computer program COUTH, "a simplified polyenergetic cell 

theory (which) is formulated to determine spatially averaged 

energy dependent thermal fluxes in the moderator, cladding, 

and fuel regions within the unit cell of a reactor lattice 

The derived spectra are then utilized in the calculations of 

the thermal integral parameters and average cross sections 

19 / 

required for reactor computations n — / Calculations at two 
temperatures provide the temperature dependence of the thermal 
parameters used here 



bi 1 ity 

The resonance escape probability temperature derivative 

can be calculated by the following scheme Weinberg and 

20 / 

Wigner- — give the- following expression for p(E), the 
resonance escape probability to energy E from source energy 



2,3,3 Resona nce Esca pe Prcba 



19/ 

— C. R, MacVean, footnote 2, p ii 

A M. Weinberg and E P Wigner, The Phys ical Theory of 
Neutron Ch ain R eactor s, p 301 University of Chicago 
Press, Chicago (1939) 



20 / 
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p(E) = exp 



-0/5) 



(E') I (E')+I, ]' i dE'/E' 



( 2 „ 30 ) 

where £ is the average lethargy increase per collision, 2 

cl 

is the energy dependent absorption cross section of the cell, 

and 2 is the scattering cross section of the cell. Equation 
s 

(2„31) is rewritten in the usual form 



p - exp 



■N (52 ) 
u vs s' 



R 1 



eff 



(2.31) 



where RI cc is the effective resonance integral, and N is 
the homogenized fuel atom density in the cell 

Logarithmic di f f erentiation with respect to temperature 
results in the following expression 

(l/p)dp/d0 = -N RI ) " 1 > 1/N )dN /d0 - (1/2. )d /d0 + 

u eft ■ s' u ' u s 7 s 



(l/Rl eff )dRl eff /do: 

( 2 . 32 » 

The atomic density of the fuel varies with temperature 
however, since fuel atoms are lost from the cell only 
through axial expansion of the fuel , we have 

(1/N ) dN /d© = -e / > (2 ,23) 

u' u u v ' 

The microscopic scattering cross section is a constant, 
therefore the scattering cross section temperature deriva- 
tive is given by 
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(i/z s )ds s /de = -p s . (2.34) 

The scatterer is water and is lost from the cell by expansion 
of water, but is gained in the cell as the horizontal cross- 
sectional area of the cell increases due to aluminum 
structural expansion. The density change in the moderator 
is therefore a combination of water and aluminum expansion, 
and 



(l/Z s )dZ s /de = -P H2 o+(2/3)f3 A1 . (2.35) 

Using Eqs . (2.33) and (2.35), Eq. (2.32) can be re- 
written as follows 

(l/p)dp/d9 - -N u RI eff (eX s )- 1 [P Hz0 - | e A1 - £ P u + 

<l/RI eff )dRI eff /d9] . 

(2.36) 

Table 2.1 gives the values of the terms in Eq. (2.36) 
for the water- to- fuel ratios investigated. 

2.3.4 Fast Fission Factor 

An estimate of the temperature effect on the fast 
fission factor can be made using an empirical relation 
for the fast fission factor for UC^-water lattices as 
given by Larrimore ^ 



21 / 

— J. A. Larrimore, footnote 14, p. 145. 
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(l/2 s )d2 s /d0 = -p g . (2.34) 

The scatterer is water and is lost from the cell by expansion 
of water, but is gained in the cell as the horizontal cross- 
sectional area of the cell increases due to aluminum 
structural expansion. The density change in the moderator 
is therefore a combination of water and aluminum expansion, 
and 

(l/2 s )c£ s /de = -Ph 2 0 +(2/3)6 A1 • (2.35) 

Using Eqs . (2,33) and (2.35), Eq. (2.32) can be re- 
written as follows 

(i/ p )dp/de = -tyu^ces,)- 1 - 4 P A1 - 4 P u + 

(l/RI eff )dRl eff /d9] . 

(2.36) 

Table 2.1 gives the values of the terms in Eq. (2.36) 
for the water- to- fuel ratios investigated. 

2.3.4 Fast Fission Factor 

An estimate of the temperature effect on the fast 

fission factor can be made using an empirical relation 

for the fast fission factor for UC^-water lattices as 

21 / 

given by Larrimore:—' 



21 / 

— J. A. Larrimore, footnote 14, p. 145. 
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e-1 = (0.1565) [1+0. 875(d„ n V m /V f )+0.288 V s /V f ] _1 (2.37) 

h 2 ° 

where e is the fast fission factor, ^ is the water density, 

Q ^ 

and V is the volume of the structural materials. Logarithmic 
differentiation, assuming V S /V^ to be small, gives 



(l/e)de/d© = 5.59(V m /V f )(e-l) 2 (e) _1 |(2p H20 /3)-3 u /3] . (2.38) 



Using these values for the 1:1 core: e = 1.03, V m /V^ = 1.012, 
9 = 20°C , p u n = 2 .05x10" 4 /°C, and p = 0.3xl0" 4 /°C, the 

r^vJ U. 

temperature effect on e as calculated by Eq. (2.38) is found 

<*- 6 

to be 0.6x10 D /°C. For all other cores e is smaller than 

2 

1.03 and the (e~l) term is dominant. The fast fission 
factor temperature dependence is therefore found to be negli- 
gible . 



2.3.5 Discussion of the Temperature Derivatives of the Nucle ar 

Parameters 

Temperature derivatives of the nuclear parameters have two 
major contributions: one which is directly related to tempera- 

ture change (such as the variation of the microscopic cross 
sections for the slow group) and one which is indirectly re- 
lated to the temperature change by the effect which this tern- 
perature change has on the thermal expansion of the materials 
in the reactor. Over the investigated temperature range, the 
temperature derivatives of the microscopic cross sections change 
by about 107 o . Over the same temperature interval, the value of 
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P>Ai and are constant, and dRI^^/d© changes by about 3/ 0 . 
Therefore, there are some terms in the expression for the 
temperature coefficient of reactivity which do not change 
appreciably over the temperature interval studies. Over 
this same temperature range, the value of q changes from 
about 0.5xl0~^7°C to about 6.0xl0~^/°C. 

The above discussion suggests an expression for a of 
the form 



a = A+Bp 



h 2 o * 



(2.37) 



Calculated values of A and B are tabulated in Table 2.3 in 
the next subsection. The experiments described in Chapter 
3 showed the predicted linear dependence although the cal- 
culated and observed values of A and B did not agree for 
all water- to-fuel ratios 



2.4 NUMERICAL CALCULATION OF THE ISOTHERMAL TEMPERATEFE 
COEFFICIENT OF REACTIVITY 

22 / 

A computer program was developed by Osias-— which 
solves the neutron balance equations, calculates the fluxes 
and adjoint fluxes, and calculates the isothermal tempera- 
ture coefficient of reactivity by Eqs . (2.9)- (2 .11) for an 

22 / 

— — D. J. Osias, private communication (August, 1964). 
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infinite cylinder reactor Calculations were made at 4, 20, 
40, 60, and 80°C. These calculations represent an approxi- 
mation to the first order perturbation theory in that the 
flux shapes calculated at 20 °C were used to compute the 
temperature coefficient at all other temperatures. All para- 
meters of the two group two region problem change with 

- 4/0 

temperature on the order of 10 / C, so that changes in the 

flux shapes with temperature are expected to contribute an 
error of less than 17 0 to the calculated values of the iso- 
thermal temperature coefficient of reactivity. 

The input numbers required and their sources are listed 
in Table 2.2. Results of the computer calculation are given 
in Table 2.3 for the five water- to- fuel ratios of the 
reactor. The calculated values of a are a linear function 



COMPUTER INPUT FOR CALCULATION OF 
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CHAPTER 3 



METHODS AND RESULTS OF THE ISOTHERMAL TEMPERATURE COEFFICIENT 

EXPERIMENT 

3.1 GENERAL EXPERIMENTAL PROCEDURE 

A typical series of runs for the temperature coefficient 
experiment was made by starting at the coldest temperature and 
progressing to higher temperatures in steps of five or ten 
degrees. Smaller temperature intervals were used in the region 
where the isothermal coefficient changed sign. 

Temperature measurements were made at three times during 
a run: at startup, at approximately half power, and at maximum 
power. The average of these temperatures was assigned to the 
run . 

All control rods were at the fully withdrawn position and 
the source was removed from the core prior to collecting data 
on the stable period of the clean core. Data were collected 
over approximately three decades increase in power level . 

After the final temperature measurement, the water was 
removed from the reactor tank to the dump storage tank, heated, 
and returned to the reactor tank for the run at the next higher 
temperature . 
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3 . 2 TEMPERATURE REGULATION AND MEASUREMENT 
3 . 2 . 1 Water Heating 

When the reactor is shut down, the moderator water is stored 
in the dump storage tank located just below the reactor tank. The 
reactor tank capacity is 2,000 gallons; that of the dump tank is 
2,200 gallons. A steam heater in the dump storage tank was used 
to heat the water at the rate of about one-half degree centigrade 
per minute. Since the system lacks automatic temperature control, 
it was difficult to obtain or maintain a given temperature to 
within one or two degrees 

3 2.2 Water Chilling 

The recirculation pump for the reactor demineralizer system 

circulated the water from the dump storage tank through a 

23 / 

Worthington water chiller— (model RWW-400) and back to the 
dump storage tank. The cooling rate was approximately one -half 
degree centigrade per hour. Because of this slow cooling rate, 
a series of experimental runs was made in steps of increasing 
temperature. The lowest obtainable temperature was found to be 
about 5°C. 

3.2,3 Temperature Measurement 

Four twenty-gauge Leeds and Northrup (#58- 1 j asbestos 
insulated, fiberglass sheathed chromel -alumel thermocouples 
placed in the core through interstitiaL holes in the top grid 



23 / 



On loan from the Department of Thermal Engineering. 
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plate were used to measure the moderator temperature. Sampling 
of the temperature took place at radii and depth as given in 
Table 3 1 



TABLE 3 o 1 



THERMOCOUPLE POSITIONS 



Water- to- 


Thermocouple 


Thermocouple 


Thermocouple Thermocoup -e 


Fuel Ratio 


1 


2 


3 


4 


LI radius-^ 


1.75-TV 


2 ,5 -IV 


8 5-11 


7 ,0-IV 


depth 


36 


18 


36 


24 


2;1 radius 


4,75-11 


5,0 111 


9 ,5-111 


2 , 25- IV 


depth 


30 


24 


12 


36 


3°1 radius 


7 ,5- IV 


2, 5- IV 


7 ,5-1 


5,0-111 


depth 


29 


25 


18 


36 


4"1 radius 


8 75-11 


c/ 


2 .5 -IV 


9,75-111 


depth 


18 




18 


18 


a/ 

— Measured in 


inches from 


the center of 


the core, Roman 


numerals 


indicate 


the grid plate sector 






— / Measured in 


inches down 


from the top 


of the upper grid 


plate 



The horizontal midplane of the fueled region of the core is 
at 30 5 inches 

c/ 

— This thermocouple was located in the reactor fill pipe for thi 
water- to-fuel ratio. 



The ice junction was located in the reactor cell next to the 
reactor, and copper multiconductor cable was used to bring the 
thermocouple voltages into the control room, A precision potent^ 
ometer (Leeds and Northrup model 8662) was used to measure the 
voltages to 1" 0 001 mV which corresponds to *J" 0,025°C, 
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3 2 4 Temperature Unifo rmity 

The time required for a step temperature difference between 
the cladding and the center of the fuel rod to decay to less than 
207c was computed to be of the order of three minutes. This 
analysis was made using the standard Fourier charts for transient 
heat transfer in cylinders. Helium fills the void between the 
fuel pellet (o.d. 0,600+0 003 inches) and the cladding 
(id O^IO^O 002 inches) The effect of this void was included 
in the calculations by assuming that the fuel pellet does not 
touch the cladding This assumption results in an overestimate 
of the transient decay time inasmuch as conductive heat transfer 
from the pellet to the cladding has been neglected. The fuel 
pellet properties used in this calculation were 



The time required to take the reactor critical is several 
times the three minute interval mentioned above About fifteen 
minutes are needed to pump the water into the reactor tank from 
the dump storage tank The reactor tank fill pump was usually 
run for several minutes after the reactor tank was full in 
order to agitate the water in the reactor tank Another eight 

minutes are required to raise the control rods. Therefore, 
normal startup procedure allowed sufficient time to establish 



fuel density 
thermal conductivity 
thermal capacity 



10 6 gm/ cm 
0 07 wat t s / cm - ° C 
0 063 cal/gm-°C 
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thermal equilibrium in the reactor tank, the maximum observed 
temperature difference between thermocouples was less than ^''C . 

33 PERIOD MEASUREMENT 

Following the establishment of a good temperature profile, 
the reactor was taken on a positive period by removal of all 
control rods. With the control rods withdrawn, the reactor is 
in a clean configuration because the fuel followers in the con- 
trol rod elements are then at the same height as the fueled 
section of the normal fuel rods 

Reactivity errors resulting from positioning errors 
should be small because the differential control rod worth 
is a minimum at the fully withdrawn position The reactivity 
effect of a simultaneous one-inch misalignment of all four 
control rods was found to be about 0 . 26 , which is within the 
statistical uncertainty of the measurement,. The specifications 
for the digital indicating system for control rod position 
call for position reproducibility to within 0 05 inch. 
Therefore, misalignment of the control rods at the fully 
withdrawn position will contribute a negligible reactivity 
error, certainly less than 0„1<|. 

Time dependent flux transients associated with control 
rod motion or source removal decay rapidly so that the actual 
reactor period is within one per cent of the stable reactor 

24 / 

period after two periods have elapsed.—- To allow transients 
24 / 

L J Templin, "Reactor Physics Constants," ANL-5800, 

2nd ed „ , p. 463. Argonne National Laboratory (1963). 
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to decay, the data recording system was not started until the 
power level had increased beyond that at criticality by about 
one decade, which is more than two periods in time. 

The suitable range of periods was from 25 to 100 seconds. 
Periods less than 25 seconds are too short to measure satis- 
factorily. Lack of temperature control of the water in the 
reactor tank prevented the use of periods much longer than 
100 seconds at temperatures much different from room tempera- 
ture because temperature drifts cause corresponding reactivity 
drifts. To avoid problems associated with very long or very 
short periods, the excess reactivity of the clean core was 
kept between about 1 $. and about 25<t by adding or removing 
peripheral fuel rods . This technique was chosen rather than 
that of compensating with control rods because with the con- 
trol rods fully withdrawn, systematic errors due to control 
rod positioning and to temperature effects on control rod 
worth were avoided. 

Two BF 2 pulse counters (Hammer-20th Century 31EB70G) 

provided the data for the period determination. The outputs 

of these counters were recorded automatically by a Hamner 

25 / 

Spectrometer, a combination scaler-timer-printer. — The 
data recording cycle was adjusted so that several data points 
were recorded during the time interval of one period. The 



25 / 



S. S. Berg, footnote 1, p. 104. 
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detectors themselves were positioned in the reactor instrument 
thimbles so that saturation did not occur for power levels used. 



3.4 DATA REDUCTION 

Period measurement raw data were analyzed by the computer 

r\ s / 

code RHOFIT,— ' a least squares fit to an exponential function. 
The computed period is then fed into a subroutine which cal- 
culates the excess reactivity by the inhour equation 

6 

Pex^eff ' J. O+V)' 1 

i=l 



where : 



^ is the neutron generation time 



(3^ is the importance weighted effective delayed neutron 

th 



yield for the i group 
is the decay constant for the precursor of the i 
group 

f3 e ff is the effective delayed neutron fraction 
T is the stable period of the reactor 
P ex /P e ff t * le excess reactivity in dollars. 



. th 



Table 3.2 lists the values of the parameters used in the inhour 
equation as well as selected values of P/P e fp and T for the 
four water- to-fuel ratios used. 



26 / 



W. E. Schilling and D. H. Bryce, 
(March, 1964) . 



private communication 
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TABLE 3.2 

REACTIVITY PARAMETERS— 



W/F 


1:1 


2:1 


3 , 1 


41 


pja) 


0.0397 


0.0317 


0,0310 


0 0321 


p 2 <*> 


0.2206 


0.2125 


0.2142 


0 2147 


p 3 (%> 


0.2023 


0.1933 


0.1942 


0 1946 


p 4 (%) 


0.3553 


0.3955 


0.3945 


0 3945 


p 5 (%> 


0.1342 


0.1229 


0.1216 


0 1207 


p 6 (%> 


0.0480 


0.0440 


0.0437 


0,0435 


P e ff ( % ) 


0.7943 


0.8630 


0.8330 


0 7955 


A-^(sec -i ) 


0.0124 


0.0124 


0,0124 


0,0124 


A2(sec ■*") 


0.0305 


0.0305 


0.0305 


0.0305 


A^ (sec - "*") 


0.1111 


0.1111 


0.1111 


0,1111 


A 4 (sec -1 ) 


0.301 


0.301 


0.301 


0,301 


Ac (sec - '*') 


i — 1 

y — ( 


1.14 


T — 1 
1 1 


1.14 


— 1 

1 

o 

CD 

CO 

> vr 

r< 


3.01 


3.01 


1 1 

o 

CO 


3.01 


A /P e f f ( sec ) 0.004041 


0.005411 


0.007155 


0,008976 


Reactivity (in $) for period T (in seconds) 


T(sec) 










20 


0.2889 


0.2799 


0.2812 


0,2819 


30 


0.2307 


0.2223 


0.2234 


0.2240 


40 


0.1937 


0.1859 


0.1868 


0 ,1874 


50 


0.1676 


0.1604 


0.1612 


0,1617 


60 


0.1480 


0.1414 


0.1421 


0.1425 


70 


0.1327 


0.1265 


0.1272 


0.1276 


80 


0.1204 


0.1146 


0.1152 


0.1156 


90 


0.1102 


0.1048 


0.1054 


0.1057 


100 


0.1017 


0.0967 


0.0971 


0,0974 



a/ 

— S. S. Berg, footnote 1, pp. 157, 210. 



232 (corrected by 



D . H . Bryce) . 
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The difference in reactivity for two runs was divided by the 
difference in temperature to obtain values of the isothermal 
coefficient a. The average temperature 9 and the corresponding 
value of the volumetric coefficient of expansion of water 6, 



h 2 o 



for the two runs were assigned to that value of a. 0 WaS 

determined by differencing a table of temperature dependent 

27/ 



Figure B.l is a plot of ( 3 „ vs 0. Tables 

fi 2 



volumes of water 
B.l - B.4 in Appendix B give the temperature, period, excess 
reactivity, and number of fuel rods for each run and for all 



water- to-fuel ratios Values of a and 6 



^ q given in the same 



tables were used in another computer code to least squares fit 
the values of a to a linear function of 0 U „ 

ii ^ v 



3.5 EXPERIMENTAL RESULTS 



3.5.1 Temperature Dependent Excess Reactiv ity 

Figures 31 through 34 show the temperature dependent 
excess reactivities of the cores investigated Uncertainties 
in the values of p are within the size of the points of the 
graphs. Nominal water- to-fuel ratios are designated X;l~2 
(read X to one dash two) where X is the water relative volume 
and the -2 indicates the core configuration All the -2 cores 
were set up in the reactor with the removable unit cell at the 

27 / 

— C. D. Hodgman (Editor-in-Chief), Handbook of Chemistry an d 
Physics , 42nd ed . , p. 2143. Chemical Rubber Publishing 
Co., Cleveland, Ohio (1960). 
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center of the lattice. As a result, all -2 cores are off center 
with respect to the reactor tank. A minimum distance of more 
than two feet separated the closest fuel rods from the grid plate 
support structure. The temperature dependent behavior of the 
reactor is not expected to be affected by the eccentricity of 
the core in the reactor tank. 



3.5.2 Isothermal Temperature Coefficient of Reactivity 

The isothermal reactivity coefficient for each of the four 
water- to-fuel ratios is shown in Figure 3.6 as a function of 
temperature and in Figure 3.7 as a function of the volumetric 
coefficient of expansion of water. Figures 3.8-3.11 give these 
values of a vs ^ in greater detail for each core. The 
straight lines in Figures 3.7-3.11 represent the equation 



a = A+B(3 



h 2 o 



(3.2) 



where A and B were determined by the method of least squares. 
Fitted values of A and B are given in Table 3.3. 



3.6 DISCUSSION OF RESULTS 

3.6.1 Dissolved Air in the Water 

Some previous temperature coefficient work was plagued 

28 29/ 

with the formation of air bubbles on the fuel rods.~— * — 

— ^ J. R. Brown, e_t al., footnote 4, p. 40. 

29 / 

— P. W. Davison, ejt al_. , footnote 5, p. 69. 
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Figure 3.1: Excess Reactivity vs Temperature for Various Fuel Rod Loadings, 
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Fijure 32. Excess Reactivity vs Temperature for Various Fuel Rod Loadings 
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Figure 3.3: Excess Reactivity vs Temperature for a Fuel Rod Loading N = 199; 
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Figure 3.4: Excess Reactivity vs Temperature for Various 

Fuel Rod Loadings, N; 4:1-2. 
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Four Water-to-Fuel Ratios. 
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Figure 3.8: Isothermal Temperature Coefficient of Reactivity vs Volumetric 

Coefficient of Expansion of Water; 2:1-2. 
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isothermal Temperature Coefficient of Reactivity vs Voiumetr 
Coefficient of Expansion of Water; 3:1-2. 
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Figure 3.10: Isothermal Temperature Coefficient of React ivi 

vs Volumetric Coefficient of Expansion of Water; 4:1-2. 
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Figure 3.11: Excess Reactivity vs Temperature; 
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FITTED 


LEAST SQUARES 


TABLE 3,3 

VALUES OF A AND B FOR THE 

a + b % 2 o 


EQUATION i = 




A(in cjt / ° C ) 


A 


Temperature 


W/F 


B( in 10 4 <f: ) 


Interval 

(°C) 


1:1 


0.29+0.02 


-0.43 5+0 ,003 


3-70 


2:1 


0 . 53+0 ,02 


-0,367+0.006 


3 - 78 


3:1 


0.82+0.02 


-0 ,268+0 . 003 


3 - 62 


4:1 


1.24+0.03 


-0,138+0.011 


7 - 43 



These small air voids were thought to be responsible for some 
unexpected results in the measurement of the temperature 
coefficient. Two precautions were used by previous investi- 
gators One was to use deaerated water to prevent the problem 
from arising. The other solution was to start at maximum 
temperature and cool the reactor. Air solubility being a 
maximum at low temperatures, by cooling the water from high 
temperatures no dissolved air would be released . 

Two precautions were taken in the present experiment to 
prevent the formation of air bubbles on the fuel rods. First, 
the fuel rods were cleaned with acetone to remove grease and 
fingerprints which might tend to collect air bubbles. Second, 
water pumping was continued after the water level reached the 
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mouth of the overflow pipe leading back to the dump storage tank 
The resulting circulation of water in the reactor tank tended to 
flush away any air bubbles (and also speed up the approach to the 
equilibrium temperature) . 

Systematic errors due to variation in dissolved air content 
of the water were estimated to be negligible In the Cornell 
Zero Power Reactor, all water heating and cooling takes place 
external to the reactor tank. The only water temperature change 
occurring in the reactor tank is a small one (1°C or less) which 
is brought about by heating or cooling of the fuel rods and the 
reactor structure by the water as it is pumped into the reactor 
tank . 

An agitator is installed in the dump storage tank to insure 
uniform properties of the water when pumped to the reactor tank 
This mixing action aerates the water; therefore it was assumed 
that the water at each temperature was a saturated solution of 
air and water . 

"Between the temperatures of 5°C and 8°C, the density of 
water saturated with air was found to be 0»000003 q gm/ml less 
than the density of air-free water." — At temperatures above 
this range, the solubility of air in water is less and therefore 

30 / 

— E. W. Washburn (Editor-in-Chief), International Critical 

Tables of Numerical Data, Physics, Chemistry and Technology . 
1st ed . , p. 26. McGraw-Hill Book Company, Inc , New York 
(1928) . 
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the density effect is presumably even smaller. Density changes 
due to thermal expansion are large by comparison (of the order 
10“ 4 /°C) . 

Under the assumption that dissolved air is added or removed 
in the dump storage tank, the change in solubility of air in the 
water results in a density change which is of no consequence for 
this experiment, 

3,6.2 Temperature Dependent Excess Reactivity 

The excess reac tivity has been determined for four water - 
to-fuel ratios over the temperature intervals given in Table 3 3 
Figures 3 1-3.4 give temperature-dependent excess reactivities 
for various fuel rod loadings with sufficient accuracy to plan 
further temperature dependence studies Four thermocouples 
with a total worth of about -2<f; were in the lattices for these 
measurements . 

For all water- to-fuel ratios, the excess reactivity curves 
have an initial positive slope. The temperature at which the 
maximum excess reactivity was observed was 8.8, 15.5, and 30. 2 
for the 1:1, 2:1, and 3:1 water- to-fuel ratios respectively. 

For the 4:1 core no maximum was reached; a curve of excess 
reactivity vs temperature for this core is shown in Figure 
3 oil. Experimental data are shown up to about 45°C; for higher 
temperatures the plotted points are from a linear extrapolation 
of a vs ( 3 . For an increase in isothermal temperature from 20 C 
to 100°C, the reactivity added would be 35<| Thus, for normal 
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room temperature operating conditions, a temperature increase to 

100°C will not, in itself, lead to prompt criticality Since 

31/ 

the Safeguards Analysis- — was based on the 4:1 core and assumed 
prompt criticality, there is no significant change in the 
analysis brought about by the determination of a positive iso- 
thermal temperature coefficient of reactivity for this core 

3.6.3 Isothermal Temperature Coefficient of Reactivity 

The isothermal temperature coefficient of reactivity has 
been measured for water- to-fuel ratios of 1:1, 2:1, 3:1, and 
4:1 over temperature intervals listed in Table 3.3. The 
accuracy of the individual measurements of a is about 0.03- 
0.06(^/ o C. The accuracy determined for a from the least squares 
fitted curve is 0.02<j:/ o C for three lattices, and 0 03<j:/ o C for 
the 4:1 lattice. 

Within the accuracy of the measurements, A and B are 
different constants for each water- to-fuel ratio Theoretical 
and experimental values of A and -B are plotted vs water- to- 
fuel ratio in Figure 3.12. The experimental curves shown 
merely connect the measured points smoothly and have no 
theoretical basis. However, they permit estimates by inter- 
polation for the 1>5:1 water- to-fuel ratio: A = 0 . 40 ^0 02 in 
4/°C and B = -0.405+0.006 in 10 4 <^ . 

31 / 

— S. S. Berg, "Final Safeguards Report to the U. S Atomic 
Energy Commission for the Cornell University Zero Power 
Reactor," CURL-4, p. 7-3. Cornell University Reactor 
Laboratory (1962) . 




for the Equation a = A+Bp„ n as a Functi- 

it 2 

Fuel Ratio. 



W'-or 
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Theory predicts and experiment verifies a linear relation- 
ship between a and however, Figure 3.12 shows disagree- 

ment between the measured and calculated values of A and B. 
This is a well-defined experiment. All parts of the reactor 
were at the same temperature. The control rods were fully 
withdrawn, and the fuel followers provided a clean core 
For each temperature, a stable reactor period was measured 
There were no large perturbations introduced by experimental 
devices. The experiment does not interfere with itself, as 
do flux measurements using foils. The conclusion therefore 
must be that the theory fails to predict the proper results, 
and further refinements in the calculations are necessary 

The only available previously published data which were 
complete enough and covered a large enough temperature range 
to permit useful plots of the isothermal temperature 
coefficient of reactivity as a function of temperature, and 
as a function of the volumetric coefficient of expansion of 

water, were for the SPERT-I reactor, a highly enriched pool 

32/ 

reactor.—— The temperature range covered was 15-97°C- 
Figure 3 13 shows both a vs 9, and a vs Fu n f° r the 

ri ^ v/ 

SPERT-I reactor. The temperature coefficient of reactivity 
is seen to be a linear function of the volumetric co- 
efficient of expansion of water for this reactor also 



54- 




Figure 3.13: Isothermal Temperature Coefficient of Reactivity vs Volumetric Coefficient 

of Expansion of Water, and vs* Temperature for the SPERT-I Reactor, 
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3 . 7 CONCLUSIONS 



The excess reactivity has been measured as a function 

of temperature for water- to-fuel ratios of 1:1, 2:1, 3:1, 

and 4:1 with sufficient accuracy to allow planning of further 

temperature dependence experiments . 

The isothermal temperature coefficient of reactivity 

has been measured for the same water- to-fuel ratios and was 

found to be a linear function of the volumetric coefficient 

of expansion of water for all lattices. Furthermore, the 

results of SPERT-I also exhibit this linearity of a vs fL, n . 

h 2 u 

Theory disagrees with measured values of A and B, indicating 
a need for more sophisticated calculations for the isothermal 
temperature coefficient of reactivity. 



CHAPTER 4 



DISADVANTAGE FACTOR MEASUREMENT THEORY 
4.1 THEORETICAL DISADVANTAGE FACTORS 



Theoretical thermal disadvantage factors are calculated 
by taking ratios of spatially averaged thermal neutron den- 
sities, which include neutrons of all energies below the 
thermal neutron cutoff energy E £ (or below the corresponding 

velocity v c ) . 

33 / 

Honeck — gives the following definitions for the thermal 

neutron density, N(r); thermal neutron density disadvantage 

factor, 6 ; and the thermal flux disadvantage factor, 6 • 
n v 

-C 



6 = 
n 



N(r ) = 



N(r, v)dv 



0 



(1/V m ) J dr N m (r) 



mod 



(1/V f ) y dr N f (r) 
fuel 



(4.1) 



(4.2) 



where 



5 = 

n 



o = 



<P 



<P 



if"/# 



9 m /? f 



HTlr /““f 

v 5 n /v 



V 



V 



V = 



0 



C f c 

vN(v)dv/ J N(v)dv 
0 



(4.3) 

(4.4) 

(4.5) 

(4.6) 



33/ 



H. C. Honeck, footnote 10, p. 50 
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in the fuel and moderator regions, and where: 




r is the position vector; 

m -P 

V , V L are the moderator and fuel volumes; 

Ny N l are the spatially averaged thermal neutron densities 
in the moderator and fuel regions; 

“HTl ~p 

cp , cp are the spatially averaged total thermal fluxes in 
the moderator and fuel regions; and 

v , v are the velocities averaged according to Eq. (4.6) 
in the moderator and fuel regions . 




thermal neutron densities, fluxes, average velocities, 
density disadvantage factors and flux disadvantage factors 
for the five water- to-fuel ratios of the Cornell University 
Zero Power Reactor. His results will be given later in 
this chapter as their use is required. 



lated to the neutron density disadvantage factor is here 
called the dysprosium, manganese, or vanadium disadvantage 
factor depending on the foils used. This quantity is defined 



4.2 EXPERIMENTAL DISADVANTAGE FACTORS 



The experimentally measured quantity § e which is re- 



by 



34 / 



C. R. MacVean, footnote 2, pp . 143-145. 
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count rate per unit volume in the moderator foil 

6 = — — — — — — . 

ex P count rate per unit volume in the fuel foil 

(4 7) 



For a foil, the count rate per unit volume C is given by 

oo 



C = e a s (1/V) K(t) 



/ d -/ 

0 



dE 2 a (E)cp(E,r) 



(4.8) 



where 



e 

a 

s 

V 

V E > 

cp(E,r) 
K ( t) 



is the counter efficiency; 
is the foil self-absorption correction; 
is the foil self-shielding correction; 
is the foil volume (integration over r_ is over the 
foil volume) ; 

is the energy dependent activation cross section; 
is the energy and position dependent flux; and 
is a time dependent function which takes into account 
the irradiation time, the time since reactor shut- 
down, and the counting time. For any given dis- 
advantage factor measurement, K(t) will be the same 
for both moderator and fuel foils and will cancel 
out of the expression for the disadvantage factor. 
In practice it need not be calculated. 
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The counter efficiency, e, will be different for the 
differently shaped foils. This effect and an experimentally 
determined correction is discussed in the next chapter. For 
purposes of this chapter, e is assumed to cancel. 

For small source- to-counter distances, the self- 
absorption factor, a, is a function of the source thickness 
35/ 

only. — Any sensitivity to geometric shape is determined 

experimentally and included in the counter efficiency 
correction factor discussed in the next chapter. For purposes 
of this chapter, then, the self -absorption factors will cancel 
in the expression for the disadvantage factor. 

The self-shielding factor, s, may vary slightly from the 
fuel to moderator foils, but this variation is neglected for 
this work and the self-shielding factors are assumed to can- 
cel in the expression for the disadvantage factor. Experi- 
mental determination of the self- shielding effect could be 
accomplished by measuring & e as a function of foil thickness 
and extrapolating to zero thickness. 

If the energy range is broken up into the thermal region 
(E varying from 0 to E c ), and the fast region (E varying from 
E c to ao ) , the following expression can be written for a cal- 
culated value 6 of the experimental quantity 5 : 

C3.1C 0 Xp 

35 / 

— W. J. Price, Nuclear Radiation Detection , p. 132. McGraw- 
Hill Book Company, Inc., New York, Toronto, and London 
(1958) . 
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oo 



e a s 



K(t)(V m ) -1 dE2 a (E)cp m (E,r) + f dE2 a (E) cp m (E , 

0 E_ 



calc 



oo 



e a s 



K(t)(V f ) _1 ^r dE2 a (E)cp f (E,r) + f dE2 fl (E) cp f (E , r>l 



0 E 

(4.9) 

On canceling e, a, s, K(t), and on using the following relation 
for the energy dependent, spatially averaged flux. 



9(E) 



= V -1 /drcp(E,r) 



(4.10) 



Equation (4.9) can be written as 



oo 



/ dE2 (E)9 m (E) + / dE2 Q (E)9 m (E) 



0 



calc 



oo 



(4.11) 



/ dES a (E)cp f (E) + f dES a (E)9 m (E) 



0 



4.2.1 1/v Thermal Cross Section Foils 

For a foil with a 1/v thermal activation cross section, 

2 (E) can be written 

2 (E) = 2 v /v (4.12) 

a v ' o o v ' 

where the subscript o indicates a reference energy, and 2 o v q 
is a constant. The calculated disadvantage factor then 
becomes 
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oo 



J dES o v o v~V(E) + f dE2 a (E)^p m (E) 



0 



calc,l/v E 



oo 



/ 'E.V'V® + / dEZ a (E)cp f (E) 



0 



(4.13) 



Using the relations cp/v = N, cp = C./E for E > E (where C- 

J c J 

is a constant), and 



v 



C 

f N(v)dv = J N(E)dE 
0 0 



(4.14) 



it follows that 



GD 



^alc , 1/v 



Z v N 1 " + C™ f dE E _1 2 (E) 
o o j J a 



oo 

2 v N f + C ? / dE E -1 2 (E) 

o o j * a ' 



(4.15) 



By simple manipulations it can be shown that 



6 i i / — & 

calc , 1/v n 



1 + cf Rl/(v a if 1 )! 

I — °-f - 

1 + C. Rl/(v o 1ST) 
j v o o 'J 



= 6 W 
n 



(4.16) 
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where 



& n is defined by Eq. (4.3); 
m f 

Cj , Cj are the epithermal flux proportionality constants; 

W is defined by the expression in the brackets; 

RI is the infinite-dilution resonance integral for 
the foil; 

v q is the thermal reference velocity (2200 m/sec); and 

a is the activation cross section of the foil at the 
o 

energy corresponding to v . 

If the second term in both the denominator and numerator 
of Eq. (4.16) were small compared to 1 , 6 ca ^ c q/ v would be 
equal to 6 . 

The resonance integral for the foil is defined by 

oo 



Rl = 



/ 



E _1 o a (E)dE 



(4.17) 



For the experimental values of RI used in the present cal- 
culations, E was 0.4 eV. However, the values of N, 6 , 
and 5^ calculated by MacVean are for a cutoff energy of 

1.6 eV. Corrections for the different values of E were 

c 

made according to the following scheme. 

Let N-^ ^ be given by 




1.6 

/ 



0 



dE N(E) o 



(4.18) 
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A similar expression defines Nq Then 

1.6 



N 0.4 - n 1.6 



/ 

0.4 



dE N (E) . 



(4.19) 



In the region 0.4 < E 1.6, let 

N(E) = cp(E)/v(E) 

= CjE“ 1 v(E )" 1 



(4.20) 



where C. is selected to be that value given by Ecp(E) from the 
36 / 

COUTH spectra—' at 0.55 eV . For v(E) in m/sec, and E in eV, 
the relationship between v and E can be written 



0.723xl0" 6 v (E) = E 1/2 



(4.21) 



Using Eqs . (4.20) and (4.21), Nq ^ is seen to be 

Nq. 4 “ »1.6- L142x10 ' 6 C j < 4 - 22 > 

2 2 

where N is in n/cm and is in n/cm -sec. 

Table 4.1 gives the values of a , v , and RI for the 
three foil materials used. 

Values of W, C j , Nq ^ , and the calculated 6 ca ^ c 
are listed in Table 4.2 for the foils used. 



4.2.2 Foils with a Non-l/v Thermal Cross Section 

For a foil whose thermal cross section deviates from 1/v, 
6 . can be computed from Eq . (4.11) by integrating the 

CalC 

energy dependent cross sections over the spectra of the fuel 
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TABLE 4.1 

RESONANCE INTEGRALS AND 2200 M/SEC CROSS SECTIONS FOR Dy 164 , 

Mn, and V 





1 64 








Dy 


Mn 


V 


RI (in barns) 


482—/ 


14.3-/ 


2 b/ 


a Q (in barns) 


2600^-Z 


13. 4^ 


4.5^-Z 



a/ 

— J. L. Crandall in "Reactor Physics Constants," 2nd ed . 

(L . J. Templin, ed.) ANL-5800, p. 164, Argonne 
National Laboratory (1963) . 

— / D. J. Hughes, Pile Neutron Research , p. 139. Addison- 
Wesley Publishing Co., Inc., Cambridge, Mass., 1953. 
c/ 

— D. J. Hughes and J. A. Harvey, "Neutron Cross Sections," 

BNL-325, pp . 7-18, Brookhaven National Laboratory 
(1955). 
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TABLE 4.2 

THERMAL PARAMETERS FOR CALCULATING DISADVANTAGE FACTORS AT 20°C 

FOR 1 / v THERMAL CROSS SECTIONS 



W/F 


m 


271 


3:1 




mod fuel 


mod fuel 


mod fuel 


Cj — (n/cm^-sec) 


.7993 .7708 


.8274 .7955 


.8371 .8028 


N-^ g — ^ (in 10 ^n/cm^) 


14.58 11.65 


28.62 21.18 


41.50 29.19 


Nq ^ — (in 10 ^n/cra^) 


13.67 10.77 


27.68 20.27 


40.54 28.27 


5 n, 0 .4 


1.269 


1.366 


1.434 


W 

Dy 


0.990 


0.992 


0.993 


5 Dy,l/v 


1.256 


1.355 


1.424 


Mn 


0.953 


0.962 


0.967 


°Mn , 1 / v 


1.209 


1.314 


1.387 


W 

w v 


0.976 


0.982 


0.986 


6 V,l/v 


1.238 


1.341 


1.414 


- C. = Ecp(E) @0.55 


eV from COUTH 


Spectra: C. R. 


MacVean, 


footnote 2, pp . 136-139. 




— ! W. E. Schilling, private communication (August, 


1964) . 


c / 

— Calculated from Eq. 


(4.22). 
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and moderator. The expression for 6 ca ^ c can be simplified 

somewhat by defining appropriate average cross sections and 

utilizing the resonance integral information above. 

The spatially averaged total thermal neutron flux for 

both the moderator and fuel regions is given by 

E c 

"cp = J cp(E) dE . (4.17) 

0 

Average thermal cross sections are defined by 



dE9(E) (4.18) 

where o’ for the moderator and fuel regions are integrated 
over the respective spectra. 

It can be shown that application of the above expressions 
gives for 5 calc 

1 + C m RI/O? 9™) 

A 2 

1 + C? RI/(of 9 f ) 

J a 

- y°>I> x (4-w) 



m 



calc 



= 6 



9 -j 

a 

a 




where X is defined by the expression in brackets, 
thermal absorber, 



For a 1/v 



Sec. 4.2 - 4.3 



67 



E 



E 



I 



c 



I 



c 



c 



a 



dE °o v o cp ^ v ~ 1/ 



dEcp(E) 



0 



0 



= C V /v 
o o 



(4.20) 



and Eq. (4.19) reduces to Eq . (4.16) as it should. 

Both manganese and vanadium thermal activation cross 
sections are very nearly l/v. The thermal activation cross 
section for dysprosium, however, does vary from l/v. The 



These dysprosium cross sections and l/v cross sections for 
manganese and vanadium were averaged over the COUTH spectra 
for water- to-fuel ratios of 1:1, 2:1, and 3:1. The values 
of these averaged cross sections, values of X, cpQ ^ (total 
thermal flux for a cutoff energy 0.4 eV), and the various 
disadvantage factors are given in Table 4.3. 

4.3 TEMPERATURE DEPENDENCE OF THE DISADVANTAGE FACTOR 
MacVean calculated the thermal neutron density and 
flux disadvantage factors according to the equations of 
Section 4.1. By calculating these thermal properties of 

— ^ R. Sher, et al . , "Low Energy Cross Sections of Dy dd \" 




Nuc . Sci. Eng. 11 , 369-376 (1961). 
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TABLE 4.3 

THERMAL PARAMETERS FOR CALCULATING DISADVANTAGE FACTORS AT 20°C 
FOR NON-l/v THERMAL CROSS SECTIONS 



W/F 


ITT 


T7T 


3:1 




mod fuel 


mod fuel 


mod fuel 


N 1>6 ~ (n/cm 3 ) 


14.58 11.65 


28.62 21.18 


41.50 29.19 


^1.6 (m/sec) 


3667 3971 


3134 3446 


2945 3265 


^1 ^ (n/cm 3 -sec) 


5.3465 4.6262 


8.9695 7.2986 


12.2218 9.5305 


Cj — ^ (n/cm 3 -sec) 


0.7993 0.7708 


0.8274 0.7955 


0.8371 0.8028 


~ ^ — (barns) 

— (barns) 
~o a v - (barns) 


1735 1610 


1937 1765 


2019 1820 


9.300 8.733 


10.22 9.452 


10.60 9.711 


3.507 3.295 


3.856 3.566 


3.999 3.664 


2 

cp Q>4 (n/cm -sec) 


4.2764 3.5943 


7.8618 6.2336 


11.1011 8.4558 


6 cp,0.4 


1.190 


1.261 


1.313 


X Dy 


0.989 


0.991 


0.993 


V 


1.269 


1.371 


1.446 


^Mn 


0.953 


0.961 


0.967 


5 Mn 


1.208 


1.310 


1.386 


X V 


0.980 


0.983 


0.987 


8 v 


1.241 


1.340 


1.414 


W. E. Schilling 
— ^ C. R. MacVean, 


, private communication (August, 1964) . 
footnote 2, p. 141. 
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the reactor at two different temperatures (20° and 40°C), the 
temperature dependent behavior of these properties was pre- 
dicted. Table 4.4 gives the calculated values of the neutron 
density disadvantage factor at 20°C, the logarithmic slope 
(l/5 n ) (dB^/dQ) , and the slope d& n /d0. 





TABLE 4.4 

TEMPERATURE DEPENDENCE OF 6 

n 




W/F 


6 

n 


< 1 /8n)( d 8 n / de >- / 

in 10" 4 /°C 


d6 n /d0 
in 10“ 4 /°C 


1:1 


1.252 


-2.76 


-3.46 


2:1 


1.351 


-4 .46 


-6.02 


3:1 


1.422 


-5.56 


-7.91 


^ c. 


R. MacVean, footnote 


2, pp . 141-143 . 





CHAPTER 5 



EXPERIMENTAL METHODS AND RESULTS FOR THERMAL DISADVANTAGE 

FACTORS 

5.1 REACTOR DESCRIPTION 
5.1.1 The Unit Cell 

The lattices of the Cornell University Zero Power 
Reactor consist of an upper fuel rod support grid plate, 
a lower fuel rod positioning grid plate, structural support 
members, the reactor tank, the fuel rods, and four control 
rod units each containing a cluster of three control ele- 
ments in an equilateral triangle of the same pitch as the 
rest of the core. These control rod clusters may be 
positioned in any location of the reactor because the 
control elements have the same outside diameter as a normal 
fuel rod. The upper portion of each control element con- 
tains the boron carbide poison. The lower portion is a fuel 
follower with the same dimensions as the normal fuel rods. 

The upper grid plate has an insert which allows the 
removal or insertion of four fuel rods as a unit in what 
is called the unit cell. These four fuel rods form a 
parallelogram of two equilateral triangles of the lattice 
pitch. For all water- to- fuel ratios, the unit cell is 
slightly off-center with respect to the upper grid plate 
and the reactor tank. All cores studied were assembled 
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with the unit cell at the center of the fuel loading to reduce 
the radial flux gradient in the unit cell. With a triangular 
(or hexagonal) lattice, one fuel rod must be the exact center 
of the core. Therefore, of the four fuel rods of the unit 
cell, one was the central element of the reactor and the other 
three were in the first hexagon out. 

5.1.2 Activation Rods 

The original design of the reactor facility included de- 
cappable fuel rods for the purpose of introducing foils into 
the fuel region of the reactor. These rods were slightly 
longer than the normal fuel rods but contained a normal 48- 
inch fuel loading. 

Because the full-length decappable fuel rod is too long 
to fit into the hood used for handling radioactive materials, 
and because it is difficult to insert and remove foils with 
this length of fuel rod, a special "activation" rod was de- 
signed. The eighteen-inch fuel section of the activation rod 
is supported by screw-in aluminum rods which make this special 
rod the same over-all length as a standard fuel element. 

(See Figure 5.1) When detached, the fuel section is short 
enough to be handled easily in the hood. It is easy to insert 
and remove foils with this rod because only a few fuel pellets 
need be handled. Because of the small size, ease in handling, 
and small number of fuel pellets handled, experiments can be 
accomplished with less radiation exposure to personnel than 
was previously possible with the longer rods. 
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Figure 5.1: 



The Activation Rod. 
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The activation rods are used for foil and wire activation 
within the fuel region of the core. They can be placed in any 
normal fuel rod position, but for the disadvantage factor work 
they were used only in the unit cell. The four fuel rod 
positions of the unit cell were filled with three activation 
rods and one normal fuel rod. The standard loading was as 
shown in Figure 5.2. 

5.2 FOIL PREPARATION 
5.2.1 Foil Types and Manufacture 

The foils used for the disadvantage factor measurement 
were five mil- thick dysprosium, manganese, and vanadium. 

The dysprosium was obtained from General Astrometals Corpora- 
tion and consisted of 5.24% dysprosium in aluminum. The 
manganese was in the form of an alloy called P-metal which 
is about 72% manganese, 18% copper, and 10% nickel. The 
P-metal was obtained from Metals and Controls, Inc., and has 
a thickness tolerance of ±0.00025 inch. The vanadium was 
obtained from H. Cross Co. 

The fuel region foils were punched from stock with a 
0 . 600±0 . 001-inch diameter die. Moderator region foils were 
punched from stock using a machined die of appropriate 
shape. Figure 5.3 shows the dimensions and shape of the 



moderator foils. 



-74- 




A ACTIVATION ROD IN CENTRAL FUEL POSITION 
OF THE UNIT CELL. 

B,C ACTIVATION RODS IN FIRST HEXAGON POSITIONS 
OF THE UNIT CELL. 

D NORMAL FUEL ELEMENT IN A FIRST HEXAGON 
POSITION OF THE UNIT CELL. 

E REMOVABLE UNIT CELL FUEL ROD SUPPORT 
PLATE. 

F NORMAL FUEL RODS IN NORMAL LATTICE POSITION. 



Figure 5.2: Top View of die Unit Cell. 
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i : i 

2: 1 
3: 1 



(in) 

0.675 

0.675 

0.675 



B 

(in) 

0.856 

1.030 
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Figure 5.3: Moderator Foil Dimensions 
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5.2.2 Foil Placement in the Reactor 

To prevent fission product contamination of the fuel 
foils, a special fuel pellet with an aluminum foil glued 
to one end with G.E. Glyptal was placed on each side of the 
fuel foils. See Figure 5.4. Ten fuel pellets were removed 
from the activation rod and an aluminum lined fuel pellet 
was inserted. The manganese fuel foil was placed next to 
the fuel pellet liner. Pellet and foil were pushed to full 
insertion with a lucite rod which was marked for trans- 
ferring the foil position to the outside of the activation 
rod. Another lined pellet was inserted, followed by a 
normal fuel pellet. This loading sequence was continued 
for dysprosium and vanadium foils until each activation 
rod had a manganese foil at the center of the rod, a 
dysprosium foil three inches below the center of the rod, 
and a vanadium foil another three inches lower. Figure 5.1 
gives the distance between the horizontal midplane of the 
fuel in the reactor and the bottom edge of the bottom fuel 
pellet in the activation rod for the three water- to-fuel 
ratios studied. 

Two moderator foils were secured to a one-eighth-inch 
thick lucite moderator support plate with double-backed 
Scotch tape. The unit cell was assembled with three acti- 
vation rods and a normal fuel rod. Using the marked lucite 
fuel foil insertion rod, the lucite moderator foil supports 
were positioned at the proper levels on the outside of the 
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Figure 5 . 4 : Fuel Foil Placement in the Activation Rod. 
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activation rods. The assembly was taped together with Band 
Aid Air Vent adhesive tape placed about three inches above 
the manganese foils and about three inches lower than the 
vanadium foils. (This tape was found to be easy to use, to 
keep its adhesive properties under water, and to come cleanly 
off the fuel rods leaving no sticky residue.) To insure 
positive positioning of the lucite moderator support plates, 
and to prevent loss of foils, waxed string was looped 
around the fuel rods immediately above and below the support 
unit. Figure 5.5 shows the assembly ready for insertion 
into the reactor 

The assembled unit cell was then loaded into the reactor 
and irradiated for a period of about ten to twenty minutes 
at flux levels corresponding to about one-half watt in power. 
Irradiation times and power levels were adjusted to give 
good count rates with minimum personnel exposure during 
handling of the fuel . 

It was necessary to remove the unit cell assembly as 
rapidly as possible to get any information at all from the 
short-lived vanadium foils (half life of 3.8 minutes). 

After a wait of about five minutes from shutdown, the acti- 
vation rods were removed to the radioactive handling hood 
for disassembly. The vanadium foils were removed, washed 
in acetone, and beta-counted on an automatic sample changer. 
About one hour after shutdown, the manganese and dysprosium 



- 79 - 



fOP EDGE OF TOP FUEL PELLET 



ADHESIVE TAPE 



9" 



tL 



FUEL FOILS, 
Mn 






a. 






FOILS 



3" 



3" 



ZzzzzZZzzS 



zzzzzzzzzz 



MODERATOR FOILS 

'A 



21 £_ 



zzzzzzzzzz 



A 



FOILS 



IZZZZZZZZZZ 



FOILS 



vv / / /-r-r^7~r 



LUCITE 
MODERATOR FOIL 
SUPPORT PLATE 



7J77/7 77 7 7 7 ' 7 



ADHESIVE TAPE 



X. 



I 



BOTTOM EDGE OF BOTTOM FUEL PELLET- 
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foils were removed from the fuel, washed in acetone, and 
counted on the automatic sample changer. 

5.3 COUNTING TECHNIQUES -- SOLID FOILS 

5.3.1 The Automatic Counter System 

38 / 

The automatic counting system described by Rushing — ' 

39 / 

and Berg — was used for the solid foil counting. There 
are 24 sample positions on a round table which pneumatically 
advances on electronic signal. The samples are held on 
pedestals which are pushed up next to a thin-window pro- 
portional flow counter for counting. The timing interval, 
table rotation, two scalers, and a printer are controlled 
electronically , 

Uniform counting geometry for the fuel foils was pro- 
vided by pedestals which were machined to accept 0.600- 
inch foils. The moderator foils were centered by eye on 
similar pedestals. 

5.3.2 Relative Fuel and Moderator Foil Counter Efficiencies 

Counter efficiency is slightly sensitive to distance 

40 / 

from the center of the pedestal. — Therefore an experi- 

38 / 

— ' J. C. Rushing, "Report of Work Done on an Automatic 
Sample Changer for Use in the Cornell University 
Nuclear Reactor Laboratory," Department of Engineering 
Physics Senior Project (1962). 

— ^ S. S. Berg, footnote 1, pp . 101-107. 
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mental geometric correction factor was determined for this 

effect as follows. Fuel and moderator foils were irradiated 

on a flux wheel in the uniform thermal flux of the TRIGA 
41/ 

hohlraum— — and both shapes were counted on the automatic 
counting system. After irradiation in a uniform flux, the 
two shapes of foils should have the same activation per unit 
volume. The foils were counted using the Geneva motion 
described in Section 5.3.4. The ratio of fuel foil count 
rate per unit volume to the moderator foil count rate per 
unit volume is the desired geometrical correction factor. 
These correction factors are given in Table 5.1 for the three 
water- to-fuel ratios studied. 



TABLE 5.1 

EXPERIMENTALLY DETERMINED COUNTER EFFICIENCY CORRECTION 

FACTORS 



W/F 1:1 


2:1 


3:1 


Factor — 0 . 990 


1.007 


1.014 


a/ 

— The statistical error is 


about 0.004 


for all factors. 



E. B. Fehr, internal memo RM-16, Cornell University 
Nuclear Reactor Laboratory (1964) . 

41 / 

— E. B. Fehr and E. P. Whitted, internal memo RM-15, 
Cornell University Reactor Laboratory (1964) . 
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5.3.3 Data Reduction for Vanadium Foils 

The short half-life of the vanadium foils made it necessary 

to analyze the data differently than the data from the other two 

foils. The vanadium foils were counted in a given order on both 

proportional counters of the automatic counting system. The 

count rate at the start of counting was about 1,000 counts per 

second. The foils were counted down to about 20 counts per 

42/ 

second. SAT I,— a local computer code, was used to correct 
each datum point back to the activity on removal from the 
reactor. These count rates on removal were then converted to 
count rate per unit volume. The moderator foil count rates 
were averaged, as were the fuel foil count rates. The ratio 
of these averaged count rates per unit volume for the moderator 
and fuel foils was corrected by the appropriate geometrical 
correction factor of Table 5.1; the result was the desired 
vanadium disadvantage factor 6^. 

5.3.4 Data Reduction for Manganese and Dysprosium Foils 
Manganese and dysprosium have half-lives of 2.59 hours 

and 2.42 hours respectively. This is sufficiently long to 
allow use of the Geneva motion counting technique. In this 
method, the foils are counted in a given order once through. 

The last foil is counted twice. The remaining foils are then 
counted in the reversed order. For a total counting time 

42/ 

— W. E. Schilling, internal memo, Cornell University Reactor 
Laboratory (1962) . 
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sufficiently short compared with the half-life of the foils, 
a simple arithmetic time average of two counts is very close 
to the actual counts at the midpoint in time. For example, 
for a total counting time of 20 minutes, the arithmetic 
average of the end point counts is within 0.3% of the true 
counts at the midpoint time. The total counting time was 
usually about 10 minutes; therefore the error is even less. 
Statistical accuracy of counting was generally of the order 
0.77, and thus the error introduced by using the Geneva 
counting method is much less than the statistical error. 

5.4 COUNTING TECHNIQUE -- DISSOLVED FOILS 
The manganese foils were counted with a "soluble" 
method in which the foils were dissolved in acid and 
aliquots of the solution were gamma-counted in a well- 
type scintillation counter. Counting of aliquots of 
solutions of foils provides two advantages over solid 
foil counting. First, the moderator and fuel foils are 
counted under identical counting geometry. Thus, there 
is no need to determine a counter efficiency correction 
factor. Second, the process of dissolving the foil 
averages the count rate of the foil over its volume. The 
count rate of the aliquot is then directly proportional 
to the total count rate of the foil itself, the constant 
of proportionality being the fraction of the total solution 
which is being counted. For both moderator and fuel foils, 
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the constant of proportionality was the same. The above total 
count rate, when divided by the volume of the foil, is exactly 
the spatially averaged count rate per unit volume needed to 
determine the manganese disadvantage factor. 

5.4.1 Foil Preparation 

After the solid foils had been counted as described in 
Section 5.3.4, they were removed to the fume hood in the 
isotope handling room. Each foil was cut with metal shears 
so that it would fit down the neck of a ten-milliter 
volumetric flask into which about five milliliters of 
3N HNO^ had previously been placed. Once the foil was dis- 
solved, the solution was diluted to 10 milliliters with dis- 
tilled water. Each solution was then mixed well. 

Three milliliters of the solution were pipetted off 
and placed into a plastic vial, and the vial was capped. 

Each foil had its own vial; each vial had three milliliters 
of solution. Two pipettes were used, one for the fuel 
foils and one for the moderator foils to prevent mixing in 
the pipette of solutions with different activities. 

5.4.2 Counting of the Solutions 

The vials fit snugly into the well of the scintillation 
counter, the top surface of the solution being below the top 
of the crystal. The geometrical arrangement in the counter 
was identical for each solution. 
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Three counting methods were initially investigated. 
First, a multichannel analyzer was used to get the gamma 
spectrum of the solutions. The area under the peak was 
determined by numerical integration. These total counts 
were divided by the foil weight to get counts per unit 
weight, and the ratio of the average counts per unit 
weight for moderator and fuel foils gave the experimental 
disadvantage factor. 

The second method was to count the solution on a 
single channel analyzer in the differential mode, adjusting 
the base line and window width to straddle the gamma peak. 
Again, the ratio of counts per unit weight for moderator 
and fuel foils was the desired disadvantage factor. 

The third method used an integral discriminator with 
the base line in the valley above the Compton edge. In 
this manner, all counts above the base line were included 
in the count rate per unit weight. 

All three methods gave the same result within experi- 
mental error. The third method was chosen because of 
simplicity of analysis, greater availability of the 
associated electronic equipment, and higher count rate 
which resulted in greater statistical accuracy. 

The aliquots were counted using the Geneva motion pre- 
viously described. Typical counting times were one minute 
per vial; typical total counts were of the order 50,000. 

The background counts were subtracted from the total counts 
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prior to averaging the counts per unit weight for the fuel 
and moderator foils. A typical data reduction is given in 
Table 5.2. 

5.5 EXPERIMENTAL RESULTS 

Disadvantage factors were measured as previously de- 
scribed at various temperatures for three water- to-fuel 
ratios. The experimental results are given in Table 5.3 

as 5 Dy’ 6 Mn’ 6 Mn,sol ( man S anese soluble), and 6 y . All 
values except the dissolved foil experiment have been 
corrected for geometry by the correction factors of Table 
5.1. Soluble manganese experiments had identical counting 
geometry and therefore needed no geometric correction 
factor . 

Values for the disadvantage factors from Table 5.3 
were least-squares fitted to a linear function of tempera- 
ture. Table 5.4 gives fitted experimental values and 
theoretical values of the disadvantage factors at 20°C. 
Both experimental and theoretical values of the slope of 
the 6 vs 9 curves are also given. 

Temperature dependent values of and are 

plotted in Figure 5.6 for the three water-to-fuel ratios 
investigated. Figure 5.7 gives 6^ vs temperature for the 
3:1-2 core. Figure 5.8 gives 6^ vs temperature for the 



same core. 
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TABLE 5.2 

DATA ANALYSIS FOR A TYPICAL MANGANESE DISADVANTAGE FACTOR 
USING DISSOLVED FOIL TECHNIQUE MEASUREMENT 

Single Channel Analyzer Settings, Integral Mode 
Fine Gain 0.5 Base Line 17.0 V 

Coarse Gain 4 Window Width 0.5V 

Clipping Time 0.5 (microseconds) Discriminator 8.5 V 



Foil—/ 

Number 


Foil 

Weight 

(gm) 


Counts 

per 

Minute 


Counts/ 

Unit 

Weight/ 

Minute 


Average 

Counts/ 

Weight, 

Moderator 

Foils 


Average 

Counts/ 

Weight, 

Fuel 

Foils 


^Mn , S o 1 


119 


0.1749 


21402 


237 719 


310 868 


238 376 


1 .304± 


307 


0.2533 


39897 


308 366 






0,007 


124 


0.1748 


21043 


239 273 








317 


0.2534 


40677 


313 370 








131 


0.1749 


21373 


238 136 









131 20629 
317 39083 
124 20774 
307 38564 
119 20527 



a/ 

— Foils numbered in the 100's are fuel foils; those in the 
300's are moderator foils. These foils are for the 
3:1 core. 

— / Counts per unit weight per minute are the sum of counts 
for each foil, corrected for background of 176 counts 
per minute, and divided by the foil weight. 
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Figure 5.6: Manganese Disadvantage Factors vs Temperature; 
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Figure 5.7: Dysprosium Disadvantage Factors vs Temperature; 3:1-2. 
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Figure 5.8: Vanadium Disadvantage Factors vs Temperature; 3:1-2. 
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5 o 6 DISCUSSION OF RESULTS 



5.6.1 Foil Construction and Positioning Errors 

Moderator foils were punched from stock using a machined 
die of appropriate shape. The radius of the circular cuts 
taken in the corners of the triangles was that of a standard 
drill size (0.6875 inch) rather than the exact radius of the 
outside of the fuel rods (0.666 inch). This gave all the 
moderator foils a slightly smaller volume than the actual 
moderator volume of interest. 

An estimate of the error associated with inexactly cut 
foils can be made as follows. First it is assumed that the 
foil thickness is uniform. A second assumption is that the 
only error in foil construction will be in making the modera- 
tor foil. This is justified somewhat in that the fuel foil 
is made from a die with less than 0,5% tolerance. The 
following equation is written for the correct experimental 
disadvantage factor for a foil with 1/v absorption cross 
section: 



where the subscript "foil" indicates a property of the foil, 
and the subscript nf indicates a property of the missing 
portion of the foil, or the "non-foil." All other symbols 



5 



correct 





(5 1) 



have been previously defined. The second assumption above 
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f f 

sets V = Equation (5 1) can then be rewritten as 



. = ° 

correct meas 



(< V foil /V "> + (^ f v: £ /^ oil v"] (5.2) 



where ^ meas is the previously reported disadvantage factor 
as given in Table 5 3 Assuming that the missing section 
of the foil is from the circular region near the fuel 
cladding, the neutron density cladding disadvantage factors 
of MacVean—' can be used to approximate This 

will result in an overestimate of the correction factor. 

Let 



R = V m /V 



m 

foil 



Then, Eq (5 2) becomes 



correc t 



meas 



R _ 1 ;i+(R-l)/o c ] 



(5,3) 



(5,4) 



where 6 is the cladding disadvantage factor 

f ITl 

Typical values of Vf Q ^ /Vf 0 ii were taken from foils 
used in the experiments Table 5 5 gives the values of 
the corrections for the three foils and three water-to- 
fuel ratios used It can be seen that the correction is 
small, and generally within the accuracy of the experiment. 
Because it is a small correction, and because it is not an 
accurately determined correction, it has not been applied 
to the experimentally measured disadvantage factors 



43 / 



C. R MacVean, footnote 2, p, 141 
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Errors resulting from improper vertical positioning of 
the moderator foils in the experimental assembly are 
expected to be negligible. The flux change in the axial 
direction at the foil position in the unit cell is of the 
order of one or two per cent per inch. Positioning accuracy 
is estimated to be about 01 inch. The resulting error is 
then seen to be about 0 ,27 c in count rate 

A small error may be introduced by the presence of 
the lucite support plate for the moderator foil. This 
plate is one-eighth- inch thick and of approximately the 
same shape as the moderator foils It is suggested that 
further experimental work be done with different thicknesses 
of lucite to investigate the magnitude of this effect. 

The error associated with the finite thickness of the 
foils has not been investigated experimentally. It was 
assumed in Chapter 4 that the self- shielding factors for 
the fuel and moderator foils were equal , In conjunction 
with the investigation of lucite thickness effect, it is 
suggested that the disadvantage factor be measured as a 
function of detector foil thickness. 

The process of dissolving the foils for the "soluble" 
disadvantage factor measurement is not expected to con- 
tribute any error to the measurement. The chemical handling 
technique employed standard volumetric flasks and volumetric 
pipettes. The pipettes are estimated to be accurate to 
about 0.25% and the flasks to about 0.1%. The error should 
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be uniform for each aliquot, and therefore cancel from the 
ratio which is the disadvantage factor. 

Weighing of the foils was accurate to 0.1 mg which 
represents a maximum error of less than 0.3% for moderator 
foils of the 1:1 core. For all other foils the error is 
less than 0.2%. 

5.6.2 Foil Induced Flux Perturbations 

There are two perturbations associated with inserting 
foils into the reactor for measurement. One effect is the 
depression of the flux by the foil itself as in the case 
of an absorber being placed into the low absorption region 
of the moderator. This would tend to reduce the measured 
value of the flux in the moderator. This effect could be 
determined by measuring the foil activity per unit weight 
as a function of foil thickness and extrapolating to zero 
foil thickness. 

The second effect is present in the fuel foil. Here, 
a highly absorbing region is replaced by a lightly absorbing 
foil. Neutrons may stream into the foil and the average 
flux in the foil would then be higher than that in a 
corresponding fuel region. This effect could also be deter- 
mined by extrapolation to zero thickness foils. 

Both of these effects tend to make the measured dis- 
advantage factor lower than the ratio of the neutron den- 
sities appropriate to the two regions under investigation. 
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An estimate of the first correction, called the 
outer flux perturbation, was made using the procedure 
outlined in ANL 5800 a correction factor, defined 

as the ratio of the activation in the perturbed flux 
to activation in the unperturbed flux was calculated 
for the three foil materials and for water- to- fuel 
ratios of 1:1 and 3:1 using the approximations given 
to the transport theory calculations of Ritchie and 
Eldridge. This correction is for a coin shaped foil, 
but the moderator foils used here were of the shape 
shown in Figure 5.3. An equivalent radius for these 
foils was taken to be 0.5 cm and 1.0 cm for the 1:1 
and 3:1 cores respectively. The correction was then 
estimated to be: 0.3-0. 8% for dysprosium, 0.6-1. 6% 
for manganese, and 0.1-0. 4% for vanadium where the 
lower values are for the 1:1 core. The correction 
increases the measured disadvantage factors by this 
amount . 

Because of the smallness of the corrections, and 
the uncertainty in their accuracy, they have not been 
applied to the measured disadvantage factors. It is 
recommended, rather, that an attempt be made to measure 
them experimentally. 



44 / 



L. J. Templin, footnote 24, p. 670. 
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5.6.3 Flux Perturbations Due to the Activation Rods 

Axial flux scans taken on the 1:1 core showed a large 
flux peaking near the ends of the fueled section of the 
activation rods. This peak is a result of replacing the 
fuel region of a normal fuel rod with the aluminum support 
pieces of the activation rod. Peaking was most pronounced 
in the 1:1 core because of undermoderation for this water- 
to-fuel ratio. The peaking may have been sufficient to 
affect the vanadium and possibly the dysprosium disadvantage 
factors in the 1:1 core. The manganese disadvantage factors 
are not expected to be influenced by this peaking because 
these foils were in the center of the fueled section of the 
activation rod. 

Peaking was less pronounced in the 2:1 and 3:1 cores. 

In addition, the fueled section of the activation rods was 
lowered in successive runs as indicated in Figure 5.1. 

The first change was made so that the vanadium foils were 
in the center of the core. This was changed again to the 
"symmetric" activation rod wherein the center of the fueled 
section of the activation rod is at the center of the fuel 
region of the core. 

The effect of the peaking may result in some inaccuracy 
in the measured value of the disadvantage factors for 
vanadium and possibly dysprosium, but the temperature de- 
pendence of the disadvantage factors should not be affected. 
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The only way to remove this effect is to resume use of the 
full length decappable fuel rods. A minimization of the 
effect can be achieved by using only one set of foils and 
placing them at the center of the fueled region of the 
activation rod. 

5.6.4 Comparison of Experiment and Theory 

Table 5.5 compares theory and experiment and shows 
that (1) dysprosium disadvantage factors lie between 4 
and 87, lower than theory, (2) manganese disadvantage 
factors lie between 2.4 and 4.5% lower than theory, and 
(3) vanadium disadvantage factors lie about 47. lower than 
theory. However, all measured values of the temperature 
dependence of the disadvantage factors agree with the 
theoretical values within the limits of experimental error 

5.7 CONCLUSIONS 

The disadvantage factor was measured for three 
water- to- fuel ratios using foils of three different 
materials. Values of the measured disadvantage factors 
are consistently a few per cent lower than the calculated 
values. Further refinements in the experimental technique, 
specifically utilization of different thickness foils and 
moderator support plates, could reasonably be expected to 
bring the experimental values of the disadvantage factor 
into closer agreement with the calculated values. 
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The temperature dependence of the disadvantage factor 
was measured for three water- to-fuel ratios. Measured 
values of the temperature dependence agree within experi- 
mental errors with the theoretical values as calculated by 
MacVean ^ 



45 / 



C. 



R. MacVean, 



footnote 2, p. 143 



CHAPTER 6 



SUMMARY AND CONCLUSIONS 

6.1 ISOTHERMAL TEMPERATURE COEFFICIENT OF REACTIVITY 

The isothermal temperature coefficient a is defined 
as dp/d© where |_ is reactivity and 0 is the temperature 
of all parts of the core. The isothermal temperature co- 
efficient was measured by differencing, at a series of 
closely spaced values of 0, values of p computed from 

45 / 

observed clean-core periods using a local code RHOFIT — 
Temperatures were determined by averaging readings from 
four thermocouples distributed through the moderator; 
auxiliary tests showed the temperature could be made 
uniform within ±0.1°C. 

A preliminary calculation of a was carried out 

using two group, two region, first order perturbation 

A 6 / 

theory as outlined by Larrimore.™ - A computer program 

calculated the fluxes and adjoint fluxes and used them 
in the standard perturbation theory reactivity equations 
to find the isothermal temperature coefficient. 

Dependence of a on temperature can be ascribed to two 
contributions i a rapidly varying one proportional to 
the volumetric coefficient of thermal expansion of water, 

45 / 

— W. E. Schilling and D. H. Bryce, private communication 
(March, 1964) . 
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and a second one, nearly constant over the investigated 
temperature range, which is due to temperature dependence 
of microscopic cross sections and other small temperature 
dependent effects. These considerations lead to an equation 
for a of the form 



a = A+Bg^ , 



( 6 . 1 ) 



where „ is the volumetric coefficient of thermal expan- 
sion for water. 

Experimental temperature coefficients were determined 
for nominal water- to-fuel ratios of 1:1, 2:1, 3:1, and 
4:1. Statistical errors in values of a are generally less 
than 107 o . For the cores studied, and over the temperature 
ranges used (1:1, 5-70°C; 2:1, 5-78°C; 3:1, 5-62°C; 

4:1, 7-45°C), the isothermal temperature coefficient is 
well represented by Eq . (6.1) where A and B are different 
constants for each core. A comparison of theoretical 
values of A and B with experimental values found from a 
least- squares fit to Eq. (6.1) shows a need for more 
sophisticated calculations to obtain better agreement with 
experiment . 



6.1.2 Conclusions 

Temperature dependent excess reactivities for various 
fuel loadings were determined for four water- to- fuel ratios 

46 / 

— J. A. Larrimore, footnote 14, pp . 185-210. 
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of the Cornell University Zero Power Reactor with sufficient 
accuracy to allow planning of further temperature dependence 
studies . 

The isothermal temperature coefficient of reactivity 
was measured for four water- to- fuel ratios with an accuracy 
which is comparable to that previously reported in the 
literature. The predicted and verified linear dependence 
of a on q had not been reported previously. 

Since the experimental results are based on period 
measurements on a clean core, and since the calculations 
were done with a simple two group model, the disagreement 
between calculated and measured values of A and B indicates 
the need for a refinement in the calculations . 



6.2 TEMPERATURE DEPENDENCE OF THERMAL DISADVANTAGE FACTORS 
6.2.1 Summary 

The neutron density disadvantage factor 5 n and the 

experimentally measured disadvantage factor & e are 

defined by the ratios 5 = N^/N^ and 6 = C 111 / (T^ where N 

J n exp 

and C are respectively the spatially averaged neutron 
density and foil count rate per unit volume; superscripts 
m and f indicate respectively the moderator and fuel 
regions of the reactor. 

47 / 

MacVean s calculation of thermal lattice parameters — 
gives the temperature dependence of the disadvantage factor. 



47 / 



C. R. MacVean, footnote 2, p. 143. 
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Since the calculated changes are small, highly accurate 
experiments were required for comparison. Experimental 
disadvantage factors were measured for three water- to- 
fuel ratios over the temperature range 5-65°C. 

An integral foil counting method was used employing 
a 0 . 600±0 .001-inch diameter foil in the fuel and an 
appropriately shaped moderator foil. Five-mil thick 
foils of dysprosium (5.24% dysprosium in aluminum), 
manganese (P-metal), and vanadium were chosen because 
their cross sections are nearly 1/v. Two counting 
methods were used. One was to dissolve the manganese foils 
in nitric ac id and gamma-count aliquots of the solutions 
in a well- type scintillation counter. The other method 
was to beta-count solid foils in a proportional flow 
counter. Results of the two methods agree within 
statistical errors. 

Experimental disadvantage factors were least-squares 
fitted to linear functions of temperature. Agreement 
between theory and experiment for the temperature depend- 
ence of the disadvantage factors is within the experimental 
error for all cores and foils, but the disadvantage factors 
themselves are a few per cent lower than the calculated 

disadvantage factors -- a result also found in the Brook- 

48 / 

haven and Westinghouse work reported by Honeck. — ' 



48 / 



H. C. Honeck, footnote 10, p. 49. 
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6.2.2 Conclusions 

Experimental determination of the thermal disadvantage 
factor was accomplished quite accurately using the technique 
of dissolving the foils and gamma-counting aliquots of the 
solutions. Beta-counting of the solid foils gave satis- 
factory results providing that measured values were corrected 
for the different geometrical counting efficiencies of 
different shapes of foils. The error for both methods was 
less than one per cent, an improvement on previously 
reported measurements. 

Theory and experiment disagree by a few per cent on 
the value of the disadvantage factor. This has been a 
common result, as mentioned above. Several factors may con- 
tribute to this disagreement: (1) flux perturbations due 

to finite thickness foils, (2) flux perturbations due to 
the lucite moderator foil support plate, (3) flux pertur- 
bations due to the use of activation rods with a shorter 
fuel region than normal fuel rods, (4) inexactly shaped 
moderator foils, and (5) a small error in the calculation 
of the disadvantage factor. The first two possibilities 
could be studied by using different thicknesses of foils 
and of moderator foil support plates. The third could be 
eliminated by using full length decappable rods, or 
minimized by placing foils only at the center of the fueled 
section of the activation rod. A more expensive die could 
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be made for punching the moderator foils to eliminate the 

fourth error source. Further study could be made on the 

calculational techniques for the disadvantage factor; 

however the estimated accuracy of 27, for the calculated 

values of the neutron density disadvantage factors by 
49 / 

MacVean — indicate that the calculations probably would 
not change much with slight refinements to the code. 

Agreement between theory and experiment for the 
temperature dependence of the disadvantage factor is 
within the experimental error; however the experimental 
error in these results varies from about 207> to about 
907, because small differences were taken between relatively 
large numbers to obtain the temperature dependence. 

6.3 PROPOSALS FOR FURTHER STUDY 
When the perturbation theory calculations of the 
temperature coefficient are brought into agreement with 
the experiments, these same calculational techniques 
can be used to predict the effect of adding soluble 
poisons to the moderator. It would then be of interest 
to measure the temperature coefficient of reactivity 
as a function of the distributed poisons. 

As indicated earlier, it would be of great interest 
to investigate the effect that changing the foil thick- 

49 / 

— C. R. MacVean, footnote 2, p. 111. 
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nesses would have on the disadvantage factor measurements. 
It might be of further interest to investigate the dis- 
advantage factor, and possibly its temperature variation, 
using resonance foils to study different energy regions. 

In conjunction with the above proposals, it would 
also be interesting to vary the thickness of the lucite 
moderator foil support plate and determine its effect on 
the measured value of the disadvantage factor. Another 
study could be made using polyethylene support plates. 

Further into the future, it might be of interest to 
measure these same parameters on cores which have different 
geometrical shapes. For example a square or long, thin 
lattice could be set up instead of a cylinder. The 
ability to properly predict the temperature coefficient 
for these geometries would be a measure of the accuracy 
of calculating the reflector contribution to the isothermal 
temperature coefficient of reactivity. 
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LATTICE PARAMETERS FOR THE CORNELL UNIVERSITY ZERO POWER 

REACTOR 

TABLE A . 1 

CORNELL UNIVERSITY ZERO POWER REACTOR FUEL ELEMENT 



PARAMETERS - 




Fuel Material 


Uranium Dioxide 


Fuel Rod Diameter, in: 


0 . 600±0 . 003 


Fuel Density, gm/cc: 


10 . 381±0 . 52% 


235 

Fuel Enrichment, Weight Percent U : 


2 .070±0 . 10 


235 

Fuel Enrichment, Atom Percent U : 


2.096+0.90 


Oxygen to Uranium Number Ratio: 


2.01:1 


Fuel Molecular Weight: 


270 . 26 (physical) 


Cladding Material: 


6061-76 aluminum 


Cladding Density, gm/cc: 


2.70 


Cladding Inner Diameter, in: 


0 . 610±0 .002 


Cladding Thickness, in: 


0 .028±0 .003 


Fuel Element Diameter, in: 


0.666+ .006 


a/ 

— S. S. Berg, see footnote 1, 
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APPENDIX B 



ISOTHERMAL TEMPERATURE COEFFICIENT DATA REDUCTION 



Tables B.l - B.4 give the following information in 
tabular form: date of the experiment, reference run 
number, number of fuel rods N, temperature 0 in °C, 
period T in seconds, excess reactivity p , reference 

numbers of the two runs used to calculate Ap/A0, Ap , A0, 
Ap/A0 = a, the average temperature for the two runs 0, 
and the volumetric coefficient of thermal expansion of 
water at the average temperature 6^ q. Figure B.l is 
a curve of Pjj q vs 0 obtained by differencing values of 
relative volume given in the H andbook of Physics and 
Chemistry ^ 



50/ 



C- D. Hodgman, 



footnote 27 , p 



2143 > 



- 111 - 



6.0 



- 112 - 




9 ( I N °C '• 

i<;ure B.l : Volumetric Coefficient of expansion of Water vs Temperature. 
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